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S, ■ ABSTRACT 
O , 

CO j In perturbation theory, the wave function of domain-wall quarks de- 

creases exponentially with the fifth coordinate. We show that, regardless 
of the quark's own momentum, the fall-off rate of the one-loop wave func- 
tion is equal to the slowest rate encountered at tree-level for any lattice 
four-momentum. We propose new domain-wall actions involving beyond- 
P-i! nearest neighbor couplings in the four physical dimensions, for which the 

r-| ! perturbative wave function decreases much faster. It is hoped that the new 

actions may preserve the good chiral properties of domain-wall fermions 
^ - up to larger values of the lattice spacing. 
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1. Introduction 



The coupling of the right-handed and the left-handed components of Wilson 
fermions through the QCD interaction leads in the continuum limit to the chiral 
anomaly [p]], but for finite lattice spacing a it also leads to lattice-artefact violations 
of chiral symmetries. This results in an additive renormalization of the quark mass, 
as well as in a severe tuning problem for four-fermion operators which are needed 
for the computation of weak matrix elements. The mass renormalization is 0(#q) 
in lattice units. Since g\ ~ (log(a Aq CD )) _1 , the mass renormalization diverges like 
(alog(a Aqcd)) -1 i n the continuum limit a — > 0. In the challenging lattice calculation 
of non-leptonic kaon decays (e.g. e'/e) the tuning problem is formidable because of 
the large number of chirality-disallowed mixings. 

In the domain- wall formulation of lattice QCD [2-6], the two chiral components 
arise as surface states on opposite boundaries of a five- dimensional lattice, and one 
expects their coupling to vanish when the size of the fifth dimension tends to infinity. 
(The fifth coordinate will be denoted s, and it takes values < s < N s — 1.) A precise 
non-perturbative characterization of chiral symmetry violations can be given in terms 
of the transfer matrix for hopping in the s-direction |3], |J . Being the result of slow 
decay of correlations in the s-direction, chiral symmetry violations are associated with 
near-unity eigenvalues of that transfer matrix. For a given gauge-field configuration 
the approach to the chiral limit is exponential iff the spectrum of the transfer matrix 
has a gap, and the fall-off rate is determined by the size of the gap. 

In full QCD there are several analytic results concerning the N s — > oo limit. The 
approach to the chiral limit is exponential in perturbation theory || [7|, ||, and the 
same is true non-perturbatively if a constrained gauge action (believed to be in the 
same universality class as the standard plaquette action) is used 0. 

For an unconstrained action one can also prove non-perturbatively that chiral 
symmetry is restored in the limit N s — > oo, provided the (finite!) number of sites 
in each of the four physical lattice dimensions is held fixed ||. When near- unity 
eigenvalues start playing a significant role, the chiral limit may be approached as 
slow as 1/N S . The proof that certain symmetries are restored in the limit N s — > oo 
is actually valid for any value of the coupling constant. But the identification of the 
restored symmetries as chiral ones depends on the fermion spectrum. It was recently 
shown that within the strong-coupling expansion the massless spectrum of the 
domain-wall lattice hamiltonian is either doubled or empty. Therefore the restored 
symmetries are not chiral at strong coupling (For further details see Appendix C.l). 

Of major importance is the question of how close to the chiral limit one gets in 
Monte-Carlo simulations. In trying to answer this question we rely on two sources. 



The first is the spectrum of the transfer matrix, or of the closely-related |3], |TT 
hermitian Wilson-Dirac operator. In the latter key finding [JT3] is that the 

spectral density of near-zero modes (corresponding to near-unity eigenvalues of the 
transfer matrix) rises by two orders of magnitude as the (quenched) coupling changes 
from Q/g 2 = (3 = 6.3 (a' 1 ~ 4 GeV) to (3 = 5.7 (a" 1 ~ 1 GeV). We hope that more 
results on the eigenvalue spectrum will be available in the future. 
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More information is available through lattice computations of various correla- 
tion functions [13-20]. A detailed numerical study of the chiral limit of domain- wall 
fermions was first carried out in the Schwinger model [I3f|. In QCD the first domain- 
wall simulations were promising, and the results for weak matrix elements (Bk, Oh) 
|14| and for the strange-quark mass [15] were in agreement with other methods. As 



of today more data is available. The pion-mass squared, which should extrapolate 
linearly to zero with the quark mass, is the most obvious measure of chiral symmetry. 
Using domain-wall fermions at quenched j3 = 6.0 (a -1 ~ 2 GeV) the extrapolated 
pion mass (for N s ~ 20) does not vanish exactly at zero quark mass but, rather, at a 
negative value of the order of few times 10~ 3 in lattice units [pl| (see also ref. ||17|| ). 
This value, however small, is in the same range as the light quark masses. 

In a sense, chiral symmetry violations are worst for the pion mass, because the 
lattice-artefact term in the PCAC relation is an ensemble average of a positive fermion 
correlator (see Appendix C.l; we comment that that lattice-artefact term is a better 
measure of chiral symmetry violation compared to the extrapolated pion mass, since 
it does not suffers from theoretical uncertainties due to chiral perturbation theory 
and due to finite- volume effects). 

No such positivity is encountered in the calculation of weak matrix elements |TT 



18], so chiral symmetry violations are expected to be smaller in this case. For example, 



in a recent simulation of four-fermion operators using the non-perturbative renormal- 



ization scheme, again at quenched (3 = 6.0, and using N s = 16, it was found |L8[ that 
mixing into wrong- chirality operators was practically zero (in comparison with 10% 
for Wilson fermions in a typical example). 

Going to a smaller value of the inverse lattice spacing, the situation at a -1 ~ 
1 GeV is unsatisfactory, as deviations from chiral symmetry are significant even for 
N s as large as 50 or 100 (in both quenched and dynamical simulations) [19|, 2(|. In 
the opposite direction, at a -1 ^ 3 GeV, no difficulties with the restoration of chiral 
symmetry have been reported [14, 12, |J. 

We believe that the existing results, especially those for weak matrix elements at 
a -1 ~ 2 GeV, do represent a breakthrough compared to the "pre domain-wall era". 
On the other hand, the results for the pion mass at a -1 ~ 2 GeV are not as good as 
one would hope for, and the present situation at a -1 ~ 1 GeV makes scaling studies 
with domain- wall fermions very difficult. 

Having summarized the situation in Monte-Carlo simulations let us return to the 
underlying physics. The key question is what mechanism(s) determine the abundance 
of near-unity eigenvalues of the transfer matrix. It is known that a few exact-unity 
eigenvalues must occur during the transition from one topological sector to another 
on a finite lattice with periodic boundary conditions [^, f3j. We believe, however, that 
the role of topology changing has been over-emphasized, for topological considerations 
alone do not explain the proliferation of near-unity eigenvalues nor the magnitude of 
the ensuing chiral symmetry violations. 

A simple explanation may be that the observed chiral symmetry violations arise 
(mainly) from generic fluctuations of the gauge field [21]. The effect of fluctuations 
need not be small, because the coupling constant used in simulations is not small 
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either. At the same time, as long as the coupling constant has not grown too much, 
perturbation theory should provide a reliable approximation of the leading quantum 
effects. 

In this paper we calculate the fifth coordinate's wave function of domain-wall 
quarks in the one-loop approximation (Sec. 2). The results lead us to consider new 
classes of domain- wall actions (Sec. 3). A preliminary account of this work was 



given in ref. | 22fl . (An alternative/complementary approach, whose relative merits 
are discussed in Sec. 4, is to employ an improved gauge action |19[| .) 

We now give an overview of the one-loop calculation. For large N s , zero bare 
quark mass and with the right-handed quark field near the s = boundary, the 
dressed fermion propagator near that boundary is 

G s ,Ap) = P+ X{s) . * v , X(s') + Reg., s, s' « N s , (1.1) 



where P± = ~(1 ±75), "Reg" stands for a continuous function of the four-momentum 

p, and T,k ~ S^-(g 2 , g 2 log(p 2 )). A unique feature of the domain-wall scheme is 

the s-coordinate wave function for (right-handed) quark modes. At tree level one has 

Xo(s)Kq s = (l-M) s . (1.2) 

(The five- dimensional mass term M is often referred to as the do main- wall height, 
and should not be confused with the quark mass 0].) By choosing M = 1 the free 
wave function can be completely localized on the boundary 

Xo (s) = hm (1 - MY = 5 S , . (1.3) 

The result of the one-loop calculation of the wave function is 

Xi(s)~s- 2 qi, (1.4) 

which contains also a power correction. Like an ordinary Wilson mass, M is renor- 
malized additively. Making an optimal choice of M we find 

g 1 = 0.5. (1.5) 

It should be noted that the difference between q and qi is 0(1). In the full one- 
loop result (eqs. ( |2.3| ) to ( |2.9| ) below) g 2 occurs as a pre-factor of relatively little 
importance. 

Let us now explain the physical origin of q±. Consider the free domain- wall 
propagator G° s / (p) for a given four- momentum p in the vicinity of the s = boundary 
at zero quark mass. The s-correlations described by this propagator are controlled 
by an exponent a(p). Each term in the propagator involves a factor exp(— da(p)) 
where d stands for either the separation \s — s'\ or the sum of distances from the 
boundary s + s'. For the standard domain- wall action one has max{exp(— a)} = 0.5 
for M — 1 where the maximum over the Brillouin zone is obtained at the "corner" 
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p n = (it, 0, 0, 0) and its three permutations. We will denote the set of global maxima 
by V. 

Now, at tree level, a fermion eigenmode with momentum p propagates indepen- 
dently of all other eigenmodes. But for any non-zero gauge coupling the fermions 
propagate in non-trivial backgrounds, and these backgrounds allow any given mo- 
mentum eigenmode to couple to all momentum eigenmodes. In particular, small- 
momentum quark modes couple to modes with 

We arrive at the following physical picture. A four- dimensional fermion mode 
created on a given s-layer mixes on that layer with the modes of V through the gauge 
field. As a mode with p 6 V, the fermion propagates with minimal suppression to 
some other layer s', where the action of the gauge field turns it back into the original 
mode. Propagation in the s-direction is therefore dominated by the modes of V 
leading to 

qi = max{exp(— ce(p))} ■ (1.6) 

For M — 1 this reduces to eq. ( |1.5| ). If we would momentarily regard the fifth direction 
as an imaginary-time direction, the above is recognized as the familiar result that 
propagation is always dominated by the lightest excitation in any given channel. The 
domain-wall case is particularly simple in that the gauge field is independent of the 
s-coordinate. 

Under certain conditions (basically that the coupling constant is not too large, 
see Appendix C.2 for a more detailed discussion) it should be possible to describe the 
results of numerical simulations too in terms of an effective wave function Xes{s) ~ 
s -i-SgS^_ rp^jg means th a t every quark's wave function is assumed to be the product 
of a four- dimensional wave function and the universal fifth-coordinate wave function 
Xeff(s). The exponential fall-off rate is accounted for by q c s. At relatively weak cou- 
pling (a -1 ^ 3 GeV) there seems to be no problem with the restoration of chiral 
symmetry, suggesting that g e ff < 1- For a^ 1 ~ 2 GeV, the rate at which chiral 
symmetry is restored depends sensitively on the observable. This, as well as other 
indications, suggest that q e g is very close to one, and the restoration of chiral sym- 
metry really follows a power-law behavior. (In ref. an estimate of q c g was given 
which, however, is unjustified because the power-law correction was ignored.) Then, 
at a -1 ~ 1 GeV the notion of a universal, localized, s-coordinate wave function breaks 
down. 

Comparing the perturbative results with the numerical data shows that, not 
surprisingly, the optimal tree-level value go = completely fails to describe that data. 
In comparison, the one-loop result qi = 0.5 lies approximately "half-way" between 
the tree-level value and the close-to-one values of q e g which seem to account for the 
results of simulations. Thus qi gives at least some indication of the actual behavior 
of the system. 

In this paper we adopt q\ as an analytic criterion for the quality of domain-wall 
actions. In Sec. 3 we consider new families of domain-wall actions involving beyond- 
nearest neighbor coupling. We compute the resulting qi, and find that values much 
smaller than 0.5 can be achieved. Finally, in Sec. 4 we discuss the relevance of our 
results to numerical simulations. 
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Some technicalities of the one-loop calculation are relegated to Appendix A. 
Higher order corrections are briefly discussed in Appendix B. An expanded discussion 
of some non-perturbative issues can be found in Appendix C. 



2. The one-loop wave function 

In this section we calculate the one-loop wave function of domain-wall quarks, rel- 
egating some of the technicalities to Appendix A. The finite- N s tree-level propagator 
and the one-loop self energy for domain-wall fermions were calculated in ref. [|1^, [7| 
(see also ref. f|). Here we will be interested in the range 1 < s < iV s , therefore we 
can use the simpler expressions for the tree-level propagator in the limit N s — > oo [|5] . 
Assuming the right-handed quark is localized near the s = boundary, the singular 
part of the tree-level propagator G Q st (p) is 

M(2- M)(l- M) s+t , . 

P +— ~u — ■ 2- 1 

if 

For M — > 1 this becomes 

P + ^4^, (2.2) 
if 

which means that the massless right-handed fermion field is fully localized on the 
boundary layer. 

The first quantum effect beyond the free theory is an additive correction, — SM, 
to the five-dimensional mass M. It arises in a mean-field approximation, or in per- 
turbation theory from tadpole diagrams. It is well known that this effect must be 
treated non-perturbatively [p3fl . We thus discard the tadpole diagrams, absorbing 
them into the tree-level action via the replacement M — > M — 5M. 

The full one-loop wave function is 

Xi(s) = \l + 6M-M\ s + 6xi(s), (2.3) 

where Sxi{s) comes from the "setting sun" diagram only. We will find Sxi{s) by 
matching the non-analytic piece of the dressed propagators with the r.h.s. of eq. ( |1 . 1| ) . 
We are interested in the behavior of Sxi{s) when M is close to its optimal mean- 
field value (see also Sec. 3). In the calculation below we thus set M = 1 + 5M in 
tadpole-improved perturbation theory. Since the tadpole-improved free propagator is 
a function of M — 5M (and not of M and 6M separately) the resulting propagator 
is identical to the ordinary propagator with M — 1. The setting-sun diagram will 
therefore be computed using the expression for the ordinary tree-level propagator for 
M = 1. 

The setting-sun diagram is depicted in Fig. |l|. Notice that we have not amputated 
the external legs |24|]. Except for s = t = (see ( j2.2Q ), the tree- level propagator is 



not singular at p — 0. To obtain a contribution to the r.h.s. of eq. (|1 . 1| ) , at least 
one of the three fermion lines must coincide with expression (|2.2|). (The kinetic self- 



energy correction E# in eq. (1.1) arises when all three lines coincide with eq. (|2.2|) , see 
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ref. JF], and will not be discussed here any further.) Assume first that the rightmost 
propagator in Fig. [I] coincides with eq. ( J2.2|) . Since the leftmost coordinate s is by 
assumption far from the boundary, we may take the limit p — > in the expressions 
for the self-energy part and for the leftmost propagator. We thus arrive at 

Sxi(s) = 9 2 C 2 Y: G°+, E+ o = g 2 C 2 E+_ 1>0 , (2.4) 

s'>0 

where C 2 is the quadratic Casimir and 

G°± = \^P± G%(p = 0) , E± = 1 tr P ± E., t (p = 0) , (2.5) 

and £ s ,t(p) is the 1PI self-energy obtained by amputating the external legs in Fig. [I} 
In the second equality of eq. (|2.4| ) we used the explicit expression for G s j far from 
the boundary (see Appendix A). 

When substituting eq. ( |2.3| ) into eq. fll.lj ) we find another term, Sxi(s'). This 
term is obtained when the leftmost propagator in Fig. [I] coincides with eq. ( |2.2| ). 
Following similar steps we now find 

6 X i(s')=g 2 C 2 Y,- sf _ 1 . (2.6) 

Thanks to a "parity" symmetry (see Appendix A) one has £+ t = £^ s . Hence eqs. (|2.4| ) 
and (213) agree. 

It remains to compute the diagonal part of the self-energy. One can write 

Ko = /_7 -0^h + (k)exp(-sa(k)). (2.7) 

(Note that the external momentum is zero.) The s-dependence enters through the 
exponential. All other factors were lumped into h + (k) (see Appendix A for more 
details). For large s, the above integral can be computed using a saddle-point ap- 
proximation. As mentioned in the introduction, the global maximum of exp(— a) 
corresponds to the lattice momentum p w and its three permutations. In the compu- 
tation we take h + = h + (p n ) outside the integral, and expand the exponent to second 




s 




s-1 






s-1 




s 



Figure 1: The 11 setting- sun" diagram. The external lines are not amputated. The fifth 
coordinate of each point, as indicated on the first row below the diagram, corresponds 
to the case where 1/f 1 (eq. ( $.!$( )) occurs on the rightmost line, and the leftmost point 



is far off the boundary. This yields eq. (2^4). The second row corresponds to eq. (2.6). 



7 



order around p n where we define k — (ir+h\\, k±). Including a factor of four to account 
for the degeneracy of the global maximum we obtain 



54 



2 73/2 s 2 2« 



7T 



(2.8) 



For the fundamental representation of SU(3) one has C 2 = 4/3. Substituting in 
eq. (|2.4j) we finally find 

r , s 2 72 2 0.788 fl\ s 

= -9* n2 73/2 (s _ 1)a2 ._ 1 - V -jr (2) , - » 1 ■ (2.9) 

Extrapolating eq. (|2.9|) to smaller values of s suggests that as soon as (or shortly after) 
we move off the boundary layer, Sxi{s) dominates over the tree- level term in eq. ( |2.3| ). 
This is true even if g 2 is small (or had the prefactor in eq. ( |2.9|) been numerically 
small). The reason is that the relative magnitude of the two terms is proportional to 
(0.5/(1 + SM — M)) s , and since 1 + 5M - M<1 this grows exponentially fast. 



3. New actions 

We have found that the s-coordinate's wave function of domain-wall quarks is 
dominated by quantum effects. The arguments of Sec. 2 show that the broadening 
of the wave-function is controlled, in the one-loop approximation, by the maximum 
of exp(— a(p)) over the Brillouin zone. This remains true for other domain- wall 
actions unless the s-couplings are drastically changed. So, if a different domain-wall 
action yields a smaller max{exp(— a)}, namely a faster fall-off of the wave function 
at the one-loop level, it is plausible that that new action also performs better non- 
perturbatively (we return to this issue in Sec. 4). 

The standard domain-wall action contains two parameters, the domain- wall height 
M and the Wilson parameter r (which is usually set equal to one). The M-dependence 
of max{exp(— a)} was investigated in ref. ||25|| . As mentioned in Sec. 2, however, the 
additive renormalization of M must be treated non-perturbatively. The optimal value 
used in simulations (M ~ 1.8) is nicely consistent with mean-field estimates. We will 
thus assume that the numerical optimization of M corresponds to setting M = 1+5M 
in tadpole-improved perturbation theory. Again, this means that we should determine 
max{exp(— a)} using the tree-level action with M — 1. As for the Wilson parameter, 
changing its value in the standard domain-wall action turns out to have little effect 
(see below). 

We will depart from the standard domain-wall action by allowing for couplings 
not only between nearest neighbors. In view of the obvious increase in computer time 
needed for the inversion of the fermion matrix, we try to be as economic as possible 
in our beyond-nearest neighbor excursion. We allow only for coupling between sites 
x and x + np, (but not e.g. for coupling between x and x + fx + v for 11 7^ v). In 
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Figure 2a: 2cosh(a:) as a junction of x = cos(pi) and £>2 = P3 = P4 = for 
the standard domain-wall action. Results are shown for four values of the Wilson 
parameter r. In all the figures of this section, the minimum seen in the plot is also 
the global minimum over the entire Brillouin zone. (At larger values of x, 2cosh(a:) 
is monotonically increasing, and eventually it diverges for x — > 1 .) 



this paper we consider explicitly n = 2 and n = 3, namely next-nearest and next- 
next-nearest couplings in the same direction. Also the modifications will be restricted 
to the four-dimensional part of the action, leaving the coupling in the fifth direction 
intact. (Note that we are interested in achieving a fast fall-off in the s-direction; any 
attempt to generate a smoother, continuum-like, behavior in the s-direction is thus 
the exact opposite of what we are aiming for.) 

The domain-wall operators considered here will have the following general form 
for zero quark mass 

Dff- = 5 3>t D + (5 s+ljt - 5 ttt )P+ + («J._i >t - 5,, t )P- , (3.1) 

with the understanding that on a finite lattice < s,t < N s — 1. The inclusion of a 
quark mass can be done in the usual way B . The four- dimensional part of the action 
is 

Dip) = i E 7m /Gv) - r W{p) + M . (3.2) 

This equation gives the tree-level operator in momentum space. The generalized 
Wilson term W(p) is a function of cos(p M ). In the kinetic term, f{p^) is an odd 
function of its argument, which we take to be sin(p M ) times a polynomial in 1 — cos(p M ). 
Later we will give explicit expressions for W(p) and f{p^)- As explained earlier we 
set M = 1 in the tree-level action, but we will use the freedom in varying the Wilson 
parameter r. Our convention is that W(p) and r are both positive. 

For any domain-wall action of the above form, the exponents a(p) are determined 
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by 

2 coshOO)) = — — g , (3.3) 

where -B(p) = 1 — M + rW(p) and a(p) > by convention. {B(p) = rW(p) for 
M = 1; for the standard domain- wall action eq. Q3.3|) reduces to eq. (|A.8p .) Lowering 
the global maximum of exp(— ac(p)) corresponds to raising the global minimum of 




Figure 2b: Same as Fig. for the action D23 with C3 = 4/3 for two values of r. 
The optimal value is r opt = 1.45, whereas r2 = 0.5 is a "reference value" defined by 
the condition r n W n = 2 at x = cos(7r) = —1 (see Sec. 4)- Notice the difference in 
vertical scales compared to Fig. \2^. 



Some insight about the features that control min{2 cosh(a)} can be obtained 
from very general considerations. One has W = for p = 0, and in all cases one 
aims for r W > 1 at p n = (ir, 0, 0, 0). As we gradually increase p\ for to n (keeping 
V2 = P3 = P4 = 0) at some value p c = (pi c , 0, 0, 0) we will have r W — 1. Were it not 
for the f(pi) term in eq. (^.3|), at p c we would obtain cosh(a) = exp(— a) = 1, namely 
no exponential suppression at all. To avoid this dangerous situation, we would like 
to have f 2 (pi) as large as possible at p\ = p\ c . 

Another danger lurks at the (fifteen non-zero) corners of the Brillouin zone. 
There, by construction, /(p M ) = 0, and so exp(— a) = (rW)^ 1 . We will therefore 
also be interested in increasing rW at the corners of the Brillouin zone. 

As a warm-up exercise let us consider the effect of varying r in the standard 
domain-wall action. In this case W = Su(l ~ cos (iV)) anc ^ f(Pn) = srn (P/J- For 
r < 1, min{2 cosh(a)} occurs at p^. We can increase rW at p v by increasing r. But 
in that case the value r W = 1 will occur at a smaller pi, where sin(pi) is smaller. 
As can be seen from Fig. ^ there is a transition region around r ~ 1.2. For larger 
values of r, min{2 cosh(a)} moves towards the point where r W = 1. As an example, 
for r = 2.0 one has r W = 1 at x = cos(pi) = 0.5, and min{2 cosh(a)} is at x ~ 0.42. 
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The largest value of min{2 cosh(a)}, obtained for r ~ 1.3 - 1.5, is around 2.8. This 
makes little improvement over the value 2.5 obtained at r = 1. 

To avoid this dead- lock we take W to be a non-linear function of cosQo^). We 
define 

W„ = X)(l-cos(p M )r. (3.4) 

i" 

W n requires couplings of sites x and x+njl. Once such coupling have been introduced 
into the generalized Wilson term, we allow them also in the kinetic term. Further 
raising of the global minimum of cosh (a) will be made possible by choosing f(p IJl ) 
that increases faster than sin(p M ), and by adjusting the Wilson parameter r. 

We now turn to the investigation of concrete actions. The minimization prob- 
lem was solved numerically. Note that the r.h.s. of eq. ( |3.3j ) can be expressed 
as a function of = cos(p M ) only. Using the invariance under permutations of 
the four components, it is enough to look for the global minimum over the range 
— 1 < Xi < %2 < x 3 < X4 < 1. (One can also study the minimization problem 
analytically. For any fi, sin(p M ) = always satisfies the extremality condition. For 
all momentum-components where sin^) 7^ one finds a coupled algebraic equation 
in x^. In all the cases we have studied, it turned out that the global minimum was 
either of the form (p min , 0, 0, 0) or else of the form (p min , p min , p min , Pmm) ■) 

We first consider an action containing next-nearest neighbors 

£ 23 = *E IMP.) ~rW 2 + M, (3.5) 

h(Pn) = sin(p M ) [1 + c 3 (l - cos(p M ))] . (3.6) 

W2 is defined in eq. ( [3.4|j . In Table 1 we show the resulting values of min{2 cosh(a)} 
and max{exp(— a)} for several values of C3. For each C3 we looked for the best value 
of r which we denote r opt . One sees that values of max{exp(— a)} much smaller than 
0.5 are feasible. A plot of 2cosh(a) for C3 = 4/3 is shown in Fig. |2E]. Notice the 
flatness of cosh(a) for — 1 < cos(pi) < at r = r opt = 1.45. 



C3 hip) ^o P t min{2cosh(a)} max{exp(-a)} 






P~ 




1.46 


2.83 


0.414 


1/3 


P 




1.14 


3.40 


0.326 


2/3 


p + 




1.19 


4.09 


0.261 


4/3 


p + 


1 3 

w 


1.45 


5.62 


0.184 


7/3 


p + p 3 


1.98 


8.06 


0.126 



Table 1: max{exp(— a)} for the action D23 at M = 1 and for various values of 
C3 (see text for the definitions). The second column gives the first two terms in the 
expansion of / 3 (p). For each C3 we show the result for r = r opt where max{exp(— a)} 
is smallest. 
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We now digress to discuss how the present work relates to the standard "im- 
provement program" (see e.g. the review |2B|). In the study of the hadron spectrum, 
only a single parameter (the bare quark mass) in the fermion action needs to be 
tuned. Once the correct continuum limit has been established, attention is focused 
on eliminating those lattice artifacts that vanish most slowly, that is, linearly with 
the lattice spacing. However, in the calculation of weak matrix elements one has 
to first establish the correct continuum limit. This is very problematic with Wilson 
or staggered fermions because, due to the loss of full chiral and/or flavor symmetry, 
many subtraction coefficients must be tuned. Controlling those subtractions by hav- 
ing good chiral and flavor properties simultaneously is thus of higher priority than 
the removal of any other lattice error. Furthermore, in the massless-quark limit 0(a) 
lattice artifacts are automatically excluded if chiral symmetry is maintained |[4j| . In 
that sense, approaching the chiral limit using domain-wall fermions encompasses the 
standard improvement program as well. 



0.25 
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-1 -0.75 -0.5 -0.25 0.25 0.5 



Figure 2c: Same as Fig. for the action D 35 with C5 = 2 for two values of r. The 
optimal value is r opt = 1.30, whereas the "reference value" is r 3 = 0.25. Notice the 
flatness of cosh(a) for —l<x< 0.25. 



Coming back to the new domain-wall action, since the new Wilson term W% starts 
off at order p 4 , the first lattice deviation from a relativistic (tree- level) dispersion 
relation comes only from the kinetic term. This is shown in the second column of 
Table 1. We observe that while increasing C3 from zero to 1/3 improves the dispersion 
relation, the opposite is true for C3 > 1/3. Although the error is formally of order a 2 , 
it might become significant if C3 is too large. To gain some idea on the magnitude 
of the error consider, say, p 2 ~ (400 MeV) 2 , which is relevant for kaon physics, on a 
lattice with a -1 ~ 2 GeV. This means a 2 p 2 ~ 1/25. For the last two rows of Table 1, 
the effect is 2% and 4% respectively. 

If next-next-nearest neighbors in the same direction are also allowed one can 
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decrease max{exp(— a)} further while maintaining a vanishing p 3 term. Let 



D a5 =i y £7 li f 5 (jp ll )-rW 3 + M. i 



(3.7) 



where again W n is defined in eq. (|3.4| ) and where 
h{Pn) = sin(p^ 



1 



1 + ^(1 - cos(p M )) + c 5 (l - cos(p M )) 2 



(3.8) 



Some values of min{2 cosh(a)} and max{exp(— a)} are shown in Table 2. A plot of 
2cosh(a) for C5 = 2 is shown in Fig. |2cj. Even for the last row in Table 2 (05 = 50), 
the deviation from Lorentz covariance is at the level of (50/4) (a 2 p 2 ) 2 ~ 2% for p 2 ~ 
(400 MeV) 2 . For c 5 = 5, the deviation is below 2% up to (700 MeV) 2 , and so on. 



c 5 


hip) 


r opt 


min{2 cosh(a)} 


max{exp(— a)} 


1 


p + §p 5 


0.87 


6.85 


0.149 


2 


p + ^p 5 


1.30 


10.50 


0.096 


3 


P+foP 5 


1.66 


13.36 


0.075 


4 


p + ip 5 


1.96 


15.74 


0.064 


5 


P + l§p 5 


2.21 


17.74 


0.057 


10 


p + 2 Mp 5 


3.18 


25.44 


0.039 


30 


P + 7j^P 5 


5.70 


45.51 


0.022 


50 


P + I2^p 5 


7.48 


59.84 


0.017 



Table 2: max{exp(— a)} for the action D35 for various values of C5 at M — 1 and 
t = r o P t- The second column gives the first two terms in the expansion of fs(p). 



We conclude with a number of technical comments. Replacing the four-dimensional 
part of the domain- wall action by D23 (-D35) approximately doubles (triples) the num- 
ber of entries in the fermion matrix. Therefore one should expect a corresponding 
increase in the cost of a single inversion of the fermion matrix at fixed N s . 

In the continuum limit, both the standard domain- wall action and the new actions 
discussed above support a single quark (one Weyl field on each boundary) for |1— M\ < 
1. In the case of the standard action there is a four-quark zone (corresponding to 
the corner p n and its permutations) for |3 — M\ < 1. When the Wilson term W n is 
employed instead, the four-quark zone is at |1 +r 2 n — M\ < 1. An additional benefit 
of the new actions is that the four-quark and the single-quark zones are separated by 
a large gap (as a function of M) for r ~ r opt . We expect that a clear gap should be 
found in simulations too, even though its precise location will likely be different from 
the weak-coupling limit. 

For the standard domain-wall action, the optimal value of M used in simulations 
agrees well with the mean- field estimate of 1 + 5M. One obtains 5M by substituting a 
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mean value u for each link variable in the Wilson term. We will assume that the new 
actions are gauged in the simplest way, namely using only products of link variables 
along straight lines. (E.g. the sites x and x + 2(1 are connected via U X ^U X+ ^ etc.) 
Using a mean link u ~ 0.8 at (3 = 6.0 (see e.g. ref. []T7]), the mean-field estimate is 
2r(3 - Au + it 2 ) ~ 0.9 r for D 2S and r(10 - 15m + 6m 2 - w 3 ) ~ 1.3 r for D 35 . 

Last, for any domain- wall operator with the form of eqs. ( [3.1|) and (|3.2| ), the 
N s — > oo limit defines an overlap-Dirac operator [27, 
relation (for a review see ref. [p8fl ) given by 



obeying the Ginsparg- Wilson 



£>gw = 1-75 e(75-D) 



(3.9) 



where e(x) = ±1 is the sign function acting on each of the eigenvalues of 75-D. 



4. Discussion 

In this paper we showed that the one-loop wave function of domain-wall fermions 
behaves like s~ 2 qf, where qi is determined by the free-fermion action. This is true 
for a wide class of domain- wall actions, including those defined in eq. fl3.1|). For the 
standard action qi = 0.5, whereas the addition of beyond-nearest neighbor couplings 
allows for much smaller values of q\. 

Only in the weak-coupling limit does q\ fully control the wave function. For 
finite g 2 up to some critical value g 2 c we expect the wave function to be proportional 
to some (up to power corrections) with q pt = q p t(g 2 )- The g 2 dependence can be 
parametrized in various ways. In Appendix B we consider the role of higher-order 
diagrams, and the parametrization 

Qpt = Ql exp(s(^ 2 Vi+9 4 V2 + - •)) > (4- 1 ) 

is found to be natural. The Taylor expansion of exp(s g 2 rji) corresponds to a family of 
1PI diagrams of all orders, where the first ^-dependent terms are two-loop diagrams. 
(Analogous statements apply to r] 2 etc.) In terms of q p t(g 2 ), one can define g 2 by the 
condition q p t{g 2 ) = 1. The existing numerical results suggest that the (quenched) 
value of Q/g 2 is very close to 6.0 for the standard domain- wall action. 

If both qi and rji were known for a given action, one could obtain a crude estimate 
of g 2 via a linear extrapolation. We have computed only the qi values, so we can only 
conjecture what trends are likely to affect g 2 . First, in the lower rows in Table 2, 
qi is extremely small. Nevertheless, if rji is large (and positive), g 2 may end up 
being approximately the same as (or, for that matter, even smaller than) for the 
standard domain-wall action. Because of gauge invariance there are vertices that 
depend linearly on C5. Since the lower rows in Table 2 come from actions with 
numerically large values of c 5 this should, indeed, lead generically to a large rji (and 
r/2 and so on). 

It is therefore safer to focus on the first few rows in Tables 1 and 2, where one 
is less prone to the above risk. The following heuristic argument suggests that, in 
that range of parameters, the new domain-wall actions may indeed be superior to the 
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standard one. When the Wilson parameter is equal to r n = 2 1 ~ n , an action containing 
the Wilson term W n (eq. (J3.4Q) gives rise to min{2 cosh(a)} = 2.5 (corresponding to 
q\ = 0.5) at p n = (tt, 0, 0, 0) for M — 1. The behavior at r = r n is therefore 
a common starting point over which we may try to improve by increasing r. As 
discussed in Sec. 3, for the standard action min{2 cosh(a)} is relatively insensitive to 
r. Its largest value (which is 2.8) is obtained around r opt ~ 1.3 - 1.5, namely r opt is 
less than 50% above r\. In comparison, for _D 2 3 at C3 = 4/3 (Fig. ^b|) the largest value 
of min{2 cosh(a;)} is achieved at r opt = 1.45 which is approximately three times r 2 . 
For D 35 at c 5 = 2 (Fig. ^cj) the best value is r opt = 1.3 which is more than five times 

The ability to reach larger values of min{2 cosh(a)} is thus correlated with an 
enhanced sensitivity to the Wilson parameter, and with a bigger ratio r opt /r n . Now, 
while 771,772,..., might in principle grow as r increases from r n to r opt , it is clear 
that as functions of the parameters of the theory their behavior will be very different 
from qi (see Appendix B). Therefore it is plausible that there exist "windows" of 
parameters where q p t{g 2 ) is controlled primarily by the decreasing qi, implying that 
the exponential suppression holds up to a larger value of g 2 . 

In Appendix C.2 we discuss how different ways of approaching the chiral limit 
are related to different forms of the spectral-density function of the (normal-ordered) 
transfer matrix. At weak coupling one expects to have a gap, namely almost all 
eigenvalues are smaller than some Ao < 1. The gap region Ao < A < 1 is not com- 
pletely devoid of eigenvalues, but their total number is drastically smaller than just 
below Ao- One also expects a big difference between the corresponding eigenfunctions. 
Those that lie outside the gap should be continuum-like modes that spread all over 
the lattice, while inside the gap the modes should be highly localized |29|. || . 

As the coupling constant increases a qualitative change takes place. Near-unity 
eigenvalues of the transfer matrix proliferate. For the standard domain-wall action 



this is a direct consequence [pi ] 1 1 of the proliferation of approximate zero modes 



of the hermitian, four-dimensional, Wilson-Dirac operator [ 2"5 , |T2| , |[|. The change 
(which seems to take place around quenched (3 = 6.0) shows the key features of the 
phenomenon known in condensed matter as localization |J0|. Due to the randomness 
of generic gauge-field configurations, in any given part of the lattice there is a finite 
probability to find a localized (approximate) zero mode. 

Viewing the (hermitian) Wilson-Dirac operator as a hamiltonian, under its action 
the fermions can hop only a single site. But with the new do main- wall actions the 
relevant hamiltonian is 75 D (see eq. ( |3.2| ) ) . Now the fermions may hop also two (or 
three) sites when the hamiltonian acts just once on a given state. It should be more 
difficult to trap the new fermions inside a small potential well, as now they have more 
ways of escaping out of it! This consideration too suggests that the critical coupling, 
where the exponential suppression is lost, may be larger for the new actions. 

The new domain-wall actions considered in this paper carry with them an obvious 
extra cost for a single inversion of the fermion matrix. One may hope to reduce chiral 
symmetry violations also by using improved gauge actions, because the latter tend 
to generate smoother configurations. If this goal is achieved, it may be a numerically 
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much cheaper way to reduce chiral symmetry violations. Using the Iwasaki action |31| 
it was found that the residual large- N s pion-mass squared (extrapolated to zero quark 
mass) drops by about a factor of two for a -1 ~ 1 GeV However, the new residual 



pion mass is still very big. Also, in thermodynamics, the Iwasaki action did not lead 
to any noticeable reduction in the residual pion mass (for a detailed discussion of 
various improvements see the first paper of ref. [19]). 



Our analysis suggests a possible explanation why the use of improved gauge 
actions has had only a limited success. If the (tree-level) domain-wall action is un- 
changed, the one-loop wave function ( |1.4| ) is still controlled by the same value of 
qi = 0.5. Only the numerical prefactor may change (cf. eq. (|2.9|) ). Hence, in this ap- 
proximation, the exponential fall-off rate is not getting any better for improved gauge 
actions. (The same reasoning applies to the use of "fat links" . One has to be careful, 
however, because this argument ignores higher-order corrections, cf. eq. ( |4.1| ).) 

In conclusion, in the one-loop approximation even a small reduction in q\ leads 
to a dramatic suppression of chiral symmetry violations for commonly used values of 
N s . If the actual quark's wave function is (even partly) correlated with the one-loop 
one, the new actions could give rise to a significantly better chiral behavior, enough 
to justify their increased simulation cost. 
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Appendix A. Some technicalities 



In this appendix we collect a few useful formulae. The inverse tree-level propa- 
gator in momentum space is 

(G%1 = i 5 s , t J2 1, sin(p M ) + W+(p) P + + W- t (p) P_ , (A.l) 

where 

W+(p) = 5 s+ht -B(p)8 s , t , (A.2) 

B(p) = 1 - M + $3(1 - cos(p M )) , (A.3) 
n 

and = W£ s {p). We consider only the case of a zero quark mass. Also, for the 

calculation of the self-energy we set M = 1 in the tree-level action (see Sec. 2). The 
tree-level propagator was computed in ref . ]5], [13], |7j . At the corners of the Brillouin 
zone (sm(p^) = 0, all p) the two chiralities decouple in eq. ( A.l ). For p — > the limit 



is singular. But at the other fifteen corners the limit is regular, leading to 

G° s , t = (W + )-l(p) P + + (W-)-l(p) P_ . (A.4) 
Specifically at p — p n and with a semi-infinite coordinate s,t — 0, 1,2,... one has 

(W + )-l = (W-); t l = - l -6{t>s) 2-' , (A.5) 

where 9{t > s) = 1 for t > s and 9(t > s) = for t < s. 

The explicit expression for the diagonal self-energy £+ 4 of eq. ( |2.5|) is given by 
eq. (36) of ref. 0. With a slight change of notation it reads (recall that the external 
momentum is zero) 



Kt = /7t^( 4 E^ 2 (^/2))^E(cos 2 (V2)(^ + G- 



)s,t 



- d 4 k 

(27T 

- sin 2 (^/2) (W-G+) S:t + ~ sin 2 (^) (G+ + G~) s ^ , (A.6) 

where 

(G ± )- 1 (j9) = J2 sin 2 b M ) + ^ ± (p)^(p) . (A.7) 

Explicit expressions for Gf t can be found in refs. ||, 0. As explained in Sec. 2 we 
are interested in S^ for s 1. In the saddle-point approximation we set the internal 
momentum on the fermion line to p v (or its permutations) in all terms, except in the 
exponential exp(— s' a) that occurs inside Gy , which is expanded to second order 
around p n using the definition 

2 cosh(a(p)) = -j^-^ . (A.8) 
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This expansion gives rise to the integrand in eq. (|2.8|) . Since sin(p^) = for all fi at 
p n , the last term in eq. ( [A. 6 ) is zero. At p n one has W ± G T = {W^)~ l . Since the 
matrices W ± and their inverses are triangular, the second term in eq. ( |A.6|) gives zero 
too. Only the first term in eq. ( |A.(j| ) contributes, and only for fi = 2,3,4, leading to 
eq. 

We also need the diagonal tree- level propagator G 0+ (eq. (|2~5| )) away from the 
boundaries at p — 0. For M — > 1 and p — > 0, the second-order operators (H^ + W / ~) s ,t 
and (W - W + ) M tend to 5 Sjt (except (W - W + ),,t at t = s = 0). As a result = 
<^s-i,t — W~ t . (For a semi-infinite s-coordinate W~ is a right-inverse of W + '.) 

Finally we observe that the interacting domain-wall action has a generalized par- 
ity symmetry. Writing the action as 

S,X,X4 Dg gl-^ ^- X4 y 4 

{U)^s>,y m (A.9) 



s,s' ;x,y;X4,y4 



one has 
where 



Ds,s';x,y;x 4 ,y 4 (U) =7 4 As',s;-£ -f,x A ,y A (U') 74 , (A. 10) 



Ug }X4] k = ^-3-k,xi\k 1 A; =1,2, 3. (A. 11) 

This looks like an ordinary parity transformation, except that we have switched the 
fifth coordinates of i/j and ip (s and s') . The above discrete symmetry implies 

£W(p,P4) = 74S t>s (-p,p 4 )74 (A. 12) 

and when sin^) = for all fi, one has trP + £ st = tr P S t)S . 



Appendix B. Beyond one loop 

In this Appendix we consider the role of higher-order corrections. Specifically 
the aim is to show how an 0(g 2 ) correction to q\ is built. The new wave-function is 
conveniently parametrized as 

5Xpt(s) oc s~ 2 q s pt , 
q s pt = (q x expQrV))* = q[ (l + g\ x s + -(g 2 r)i) 2 s 2 + ■ ■ ■) . (B.l) 

The last expression suggests that, in the wave function Sx P t{s), the 0(g 2 ) correction 
to qi arises from a resummation of perturbation theory. If true, at any finite order 
we should find contributions to the wave function whose structure is Sxi{s) times an 
increasing power of s. (We likewise expect additional 0(g 4 ) etc. corrections to q x , cf. 
eq. ( |4.1| ); the arguments below are, however, too crude to tell how the power-law part 
of the wave function depends on g 2 .) 
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s-1 



Figure 3a: A two-loop diagram that (likely) contributes to eq. ( \B. 1 ). In this figure 
and the following ones, the fifth coordinates correspond to the case where 1/p 1 occurs 
on the rightmost line, cf. eqs. ( \2.Q ) and $T 



We believe that a term g r\ s d~Xi( s ) arises from the two- loop diagrams in Figs. |3a 
and |3b] (and not from Figs. |3c] or [53, see below). Following Sec. 2, we will consider 
for defmiteness the case where the rightmost fermion line in Fig. 3a corresponds 
to the singular part of the tree- level propagator, eq. ( |2.2j ), and therefore the fifth 
coordinate of the rightmost vertex is zero. In Fig. |3a] the momenta on all three 
internal fermion lines can be simultaneously equal (or close) to p n (recall that the 
external momentum is (close to) zero). Then, the one-loop exponent qi = \ comes 
with a power \s — t'\ + \t' — t\ + 1 which, for given s, is minimal when the points 
are ordered: s > t' > t > 0. (Here we have ignored the difference between s and 
s — 1 which is negligible for s ^> 1.) When the points are not ordered we obtain 
an exponentially convergent series in the excessive length of the fermion's trajectory. 
For simplicity we will assume that the points are ordered, as we are only interested 
here in arguing that a contribution proportional to g 2 s5xi{s) exists. (However, the 
full series has to be summed in order to obtain the numerical value of rj.) 




Figure 3b: Another two-loop diagram that contributes to eq. ( \B. 1\ ) 



Next consider the integration over the loop momentum of the inner loop. As 
explained in Sec. 2, this integration should give rise to a factor of \t' — t\~ 2 . Since 
the sum J2m m ~ 2 is convergent, as a crude approximation one can say that the points 
t and t' are forced to be close together. Once this extra constraint has been taken 
into account, the remaining expression is independent of t (and t'). Therefore the 
t-summation (approximated by an integral) gives J S dt = s. The gaussian integration 
over the momentum of the outer loop gives rise roughly to the factor s~ 2 (associated 
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with 5xi(s)) as before. We conclude that a term proportional to g 2 s5xi(s) may 
indeed arise from Fig. |3a]. (A technical complication is that, if the two loop-momenta 
are both equal to p n , the momentum flowing through the inner gauge-boson line is 
zero. The singularity (k—p)~ 2 is integrable in four dimensions, but its existence makes 
the actual calculation quite complicated. We believe that the above considerations 



are robust enough to grasp the dominant behavior of Fig. 3a 








Figure 3c: This two-loop diagram does not contribute to eq. (B.l) 





s-1 



Figure 3d: A reducible two-loop diagram that does not contribute to eq. (\B.1\) 



The above argument easily generalizes to the higher-loop diagrams of Fig. [4a[ At 
the nth-order one has n — 1 "inner" loops. The fifth coordinates are pair-wise close, 
but otherwise are constrained only by ordering. Hence one expects a contribution 
proportional to 

ps rti rt n -2 S n ^ 

dt x dt 2 ... dt n _ x = — , (B.2) 

Jo Jo Jo [n — 1)! 

in agreement with eq. (p.l| ). A similar reasoning implies that the two- loop and n- 
loop diagrams of Figs. [J5] and ^Fj respectively also contribute to the Taylor series in 
eq. flBJ. 

We now want to explain why the diagrams of Figs, [k] and [3d| do not contribute to 
the r.h.s. of eq. ( |B.1| ). Consider first the reducible diagram Fig. |3cl|. The momentum 
on the middle fermion line is zero (being equal to the external momentum). Hence 
the points t and t' are very close (see Sec. 2 and Appendix A). Ignoring the difference 
between t and t' the two gaussian integrations give rise to the product \s — t\~ 2 t~ 2 . 
The t-summation is then dominated by t-values which are either close to zero or to 
s (being the fifth coordinates of the rightmost and leftmost vertices). Therefore the 
result behaves like s~ 2 (and not like s~ 2 s = s -1 as in the case of Fig. 
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/ rx x 



s s-1 t' n , t M t; t 2 t; t, o o 

Figure 4a: ^4n n-loop diagram that (likely) contributes to eq. I \B.1\ ). 




While the diagram in Fig. 3c is irreducible, the momenta on the internal fermion 
lines cannot be all equal to p n . For example, if the momenta on the first and third 
internal lines are equal to p^, then the momentum on middle line is zero. Hence, 
again, the points t and t' must be very close together, as well as close to either zero 
or s, and the final result is proportional to s~ 2 as in the case of Fig. Bd. 



Appendix C. Some non-perturbative observations 

C.l. Chiral-symmetry restoration 

We outline here the proof of chiral-symmetry restoration given in ref. in the 
light of later works (in particular ref. |TIJ). For definiteness we focus on the 
anomalous term in the lattice PCAC relation. (Its vanishing implies that the pion 
mass will be zero after taking the infinite-volume and massless-quark limits in that 
order.) The PCAC relation using domain-wall fermions reads |J 

^ (A^x) J b 5 (y)) = 2% J b M) + 2 (J% q {x) J*(y)) + contact term. (C.l) 

Here A M is the backward lattice derivative and m the bare quark mass. A%Jx) 
is the Noether current of a lattice transformation that assigns opposite charges to 
fermions on the half-spaces < s < N s /2 and N s /2 < s < N s . This reduces to 
a chiral transformation on the quark states, as long as they are localized in their 
respective half-spaces. The pseudo-scalar density Jg(x) is composed of fermion vari- 
ables situated on the two boundaries, and serves as the standard interpolating field 
for pions. The anomalous (lattice-artefact) term in this relation involves another 
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pseudo-scalar density, J^ q (x), which is localized on two s- layers exactly half-way in 
the fifth direction. (We are assuming an undoubled quark spectrum (one quark field 
for each five-dimensional fermion field, in most applications) which is true at weak 
coupling. As mentioned in the introduction, the massless spectrum changes at strong 
coupling in which case the above transformations are no longer chiral.) 



Before we proceed with the discussion of the anomalous term we have to address 
a technical point. Let T = J2i \i) U (i\ be the spectral decomposition of the (positive) 
first-quantized transfer matrix. We define a "normal-ordered" transfer matrix Q via 
its spectral decomposition Q — J2i K) \ (*| where Aj = min {tj, t^ 1 }. Physically, the 
operation of replacing £ 4 > 1 by its inverse amounts to filling the Dirac sea. The 
spectrum of Q lies in the interval < A» < 1. One can show [TT] that the exact 



domain-wall propagator in a given background field is a sum of terms, each of which 
involves the matrix Q raised to a positive power which is a function of the fifth 
coordinates. 

Let us now assume that, for a given background field, the spectrum of Q lies 
in the interval < Aj < Ao where Ao < 1. Coming back to the anomalous term in 
eq. ( P.l|) , for a non-singlet current it involves the propagation of two fermions over 
an s-separation equal to N s /2 (or N s /2 ± 1). The anomalous term is thus bounded 
by (A^ 2 ) 2 = A^ s times a constant. (Using the second-quantized transfer matrix 
formalism one can show that the proportionality constant is finite, being the norm of 
a product of bounded operators ||.) Moreover, if Q has no eigenvalues larger than A 
for all gauge fields (a condition which is satisfied for a constrained gauge action ||) 
one finds that the anomalous term falls exponentially after the functional averaging 
over the gauge field. 

In ref. || we showed that a very weak bound on the anomalous correlator exists 
even if there is no gap at all. The point is that exact-unity eigenvalues of the transfer 
matrix (hence of its normal-ordered version, Q, too) exist only on a submanifold 
of the lattice gauge-field space defined by the condition det Dw = where Dw is 
the hermitian four-dimensional Wilson-Dirac operator. As explained above, when a 
fermion propagates across an s-separation N s /2, the propagator is bounded by (and, 
generically, falls like) A^ 2 where A < 1 is the largest eigenvalue of Q. But A^ 2 
is negligible unless Ao = 1 — 0(1/ N s ). This condition, in turn, will be satisfied only 
for gauge-field configurations whose distance from the above submanifold does not 
exceed 0(1/N S ). The volume of the (compact) gauge-field subspace contributing to 
the anomalous correlator is therefore finite, and shrinks like 1/N S , implying a similar 
bound on the correlator itself. We comment that the restoration of chiral symmetry 
is consistent with the fact that the overlap operator defined by the N s — > oo limit 
(see eq. ( |3.9| )) admits Liischer's chiral symmetry |32| (whose generators are functions 
of the gauge field). 

C.2. Spectral density and effective wave function 

In the infinite- volume limit one can define a spectral function Pq(A) associated 
with the normal-ordered transfer matrix Q introduced in the previous subsection, 



22 



whose support is (contained in) the interval [0,1]. Here we will not attempt to 
compute any spectral function. Instead, we adopt a "phenomenological" point of 
view. We will assume that a single, continuous, spectral density function p(A) has 
been defined as a suitable configuration average of Pq(A). Using the considerations 
of the previous subsection, the aim is to see how different forms of the spectral 
function lead to different ways of approaching the chiral limit. One can envisage 
three prototype scenarios which are listed below. (A recent treatment of domain-wall 



fermions based on spectral integrals can be found in ref. |33fl . A numerical study of 
a spectral quantity which is closely related to p(l) can be found in ref. |12|]. See also 
ref. 0.) 

1) Exponential suppression. Assume that the support of p(A) is the interval < 
A < A where A < 1 and that, close to A , p(A) vanishes like (A — A) 5 with 6 > 0. 
Consider the propagation of a single fermion from the boundary layer s' — to some 
other layer s (We assume 1 <C s ^ N s /2). This involves the integral 

/ A ° dX p(A) X s ~ Aq / A ° dX p(A) exp is ) ~ s- 1 - 5 X s . (C.2) 

Jo Jo \ A / 

In order to obtain the power-law correction we have used the assumed behavior of 
p(A) close to Ao, and wrote (A/Ao) s = exp(s log(l + (A — Ao)/Ao)). (The one-loop 
result of Sec. 2 corresponds to A = 1/2 and 5 = 1.) If n fermions propagate across 
a similar s-interval, we will obtain the factor s~ x ~ 5 Xq for each of them. We may 
therefore consider Xefr(s) = s~ 1 ~ s Xq as the effective s-coordinate wave function for 
all quark states. Since the anomalous divergence J$ q (x) is a fermion bilinear, chiral 
symmetry violations should fall like x^ s (N s /2) ~ N s ~ 2lyl+ ^ Xq s . (The effective wave 
function describing the propagation of n fermions could be somewhat different from 
X™ff(<s) due to interactions between the different particles. Since chiral symmetry 
violations are related to J| q (x), the relevant effective wave function is always the one 
extracted from the sector with one fermion and one antifermion.) 

2) Power-law suppression. Assume that Ao = 1 but with a vanishing p(l), namely 
p(A) ~ (1 - A) 5 for A ~ 1 with 5 > 0. In that case the result of the spectral 
integral (|C.2|) will be s _1_<5 . We might still speak of an effective wave function 
Xcs(s) = s^ 1 ^ 5 and expect chiral symmetry violations to fall like N s ~ 2( - 1+5 \ 

3) (Almost) no suppression. Last assume that Ao = 1 and that p(l) is non-zero. 
Remember now the submanifold discussed in the previous subsection of gauge fields 
supporting an eigenvalue one of Q. If p(l) is finite, configurations close to that sub- 
manifold must have a non-negligible Boltzmann weight. As we have explained, in 
this case the only suppression of long-range s-correlations comes from phase space 
considerations (the need to pick a configuration located 0(1/ N s ) away from that sub- 
manifold). As a result, chiral symmetry violations fall roughly like p(l)/N s , and the 
concept of a localized, effective wave function breaks down, (n-fermion correlations 
fall like 1/N S too, and not like 1/N s n .) 

For clarity, we have presented above the three mathematically distinct scenarios. 
In reality, however, one is likely to encounter a more complicated behavior, char- 
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acterized by the existence of a crossover region. Let us reexamine the exponential 
suppression scenario. As discussed in the previous subsection, in the ensemble of all 
gauge-field configurations (as opposed to the case where the plaquette is constrained 
to be everywhere small H) there is a non-zero probability of finding an eigenvalue 
arbitrarily close to one, for any g > 0. Consider the integrated spectral density 
X = fl Q dXp(X), and suppose that X is comparable to for some N (up to power 
corrections). For N s ^ iVo, chiral symmetry violations will fall like N s ~ 2 ^ 1+5 ^ Xq" . 
The point is that for any Ao significantly smaller than one, and iVo of the order of 
(few times) ten, will be so small that one will never have to use N s ^ iVo in a 
simulation. For practical purposes this scenario is therefore the same as the purely- 
exponential suppression scenario. If, nevertheless, very large values of N s will be 
tried, then around N s ~ No a crossover to some slower fall-off rate will be encoun- 
tered. (For instance, a crossover to a slower exponential fall-off rate has been observed 
in the Schwinger model P~5fl.) 

Last consider the relation between the domain-wall actions discussed in Sec. 3 and 
the associated overlap operators (see eq. (|3.9|)). The spectral function that controls 
the approach to the chiral limit of domain-wall fermions also controls the localization 
range of the overlap operator. As is clear from the above discussion, problems start 
when there is no gap, namely when it is not possible to identify a range Ao < A < 1 
(with Ao < 1) where the eigenvalue density is (practically) zero. In that case both the 
overlap operator and the generators of the associated Liischer symmetries become 
non-local. This is the counter-part of the loss of exponential suppression in the 
domain-wall case. 
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